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ABSTRACT 

We present a Chandra X-ray Observatory investigation of the millisecond pulsars (MSPs) in the 
globular cluster M28 (NGC 6626). In what is one of the deepest X-ray observations of a globular 
cluster, we firmly detect seven and possibly detect two of the twelve known M28 pulsars. With the 
exception of PSRs B1821-24 and J1824-2452II, the detected pulsars have relatively soft spectra, with 
X-ray luminosities lO'^""-^^ ergs s'-^ (0.3-8 keV), similar to most "recycled" pulsars in 47 Tucanae and 
the field of the Galaxy, implying thermal emission from the pulsar magnetic polar caps. We present the 
most detailed X-ray spectrum to date of the energetic PSR B1821-24. It is well described by a purely 
non-thermal spectrum with spectral photon index F = 1.23 and luminosity 1.4 x 10^'^O(D/5.5 kpc)^ 
ergs (0.3-8 keV), where G is the fraction of the sky covered by the X-ray emission beam(s). We 
find no evidence for the previously reported line emission feature around 3.3 keV, most likely as a 
consequence of improvements in instrument calibration. The X-ray spectrum and pulse profile of PSR 
B 182 1-24 suggest that the bulk of unpulsed emission from this pulsar is not of thermal origin, and 
is likely due to low-level non-thermal magnetospheric radiation, an unresolved pulsar wind nebula, 
and/or small-angle scattering of the pulsed X-rays by interstellar dust grains. The peculiar binary PSR 
J1824-2452H shows a relatively hard X-ray spectrum and possible variability at the binary period, 
indicative of an intrabinary shock formed by interaction between the relativistic pulsar wind and 
matter from its non-degenerate companion star. 

Subject headings: globular clusters: general — globular clusters: individual (NGC 6626) — pulsars: 
general — pulsars: individual (PSR B1821-24, PSR J1824-2452H) — stars: neu- 
tron — X-rays: stars 



1. INTRODUCTION 

Rotation-powered millisecond pulsars (MSPs) are faint 
X-ray sources, with typical luminosities of Lx ~ iQ^o-ai 
[ ergs s~-^ and up to ^ lO'^^ ergs s~^. Due to the long 
exposures required to study MSPs in detail, globular 
clusters provide a particularly efficient way to observe 
I them in large numbers. This has only become feasible 
in the past decade, owing to the sub-arcsecond angu- 
lar resolution of the Chandra X-Ray Observatory that 
] is essential for individually resolving the numerous X- 
ray sources in the dense cores of clusters. Thus far, 
a total of ~30 X-ray counterparts of MSPs have been 
ident ified with Chan dra in globul ar clusters, inclu ding 47 
' Tuc (|Grindlav et all [2002: Bogd anov etaLl [20061) . NGC 
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6397 (IGrindlav et all 120 021: IBogdanov et al.| [20Tol). M28 



("Bec ker et al.ll2003D. M 4 (iBassa et al.ll2004lj . NGC 6752 
(iD'Amico et al.l mm. M71 ([Eisner et all [2008[) . and 



Terzan 5 ([Heinke et al.ll2006[ ). This appreciable sample 
has already provided important insight into the general 
X-ray properties of the MSP population in the Galaxy. 

The globular cluster Messier 28 (NGC 6626) is at a 
distance oi D fa 5 .5 kpc. The reddening towards M28 
E{B ~V) ^ 0.43 ([Harrisl 119961 ) implies a hydrogen col- 
umn density of A'h ~ 2.4 x 10^^ cm~^. For over two 
decades, M28 was known to host a single milli second pul- 
sar (MSP), PSR B1821-24 (|Lvne et al.l ll987). Recently, 
the Green Bank Telescope d etected 11 additional radio 
pulsars in M28 ([Beginl [20060 making it the cluster with 
the third largest number of known pulsars to date after 
Terzan 5 (at leas t 33, see [R ansom ct al. 20 0j|) and 47 
Tucanae (23, see iCamilo et al.ll200a IFreire etall [200l 
and references thereinf^. These new pulsar discoveries 
make M28 a well-suited globular cluster for effective X- 
r ay studies of a sizab le population of MSPs. 

[Becker et al.l ([2003[ ) have previously reported on a set 
of ^^40. 9 ks Chandra X-ray Observatory ACIS-S obser- 
vations of M28 (ObsIDs 2683, 2684, 2685) which probed 
down to a luminosity of Lx ~ 6 x lO'^'^ ergs s~^ (as- 
suming a source detection threshold of ~10 counts) and 
detected 46 sources, only 12 of which are within one core 
radius (0.24') and 26 within the half-mass radius (1.56'). 

Fo r an up-to-date list of all known globular cluster MSPs see 
|http : //ww w . nai c . edu/-pf reire/GCpsr .html 
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TABLE 1 
Chandra Observations of M28 



Telescope/ 


Epoch of 


Observation 


Exposure 


Instrument 


Observation 


ID 


Time (ks) 


Chandra/ACIS-S 


2002 Jul 4 


2684 


12.9 


Chandra/ ACIS-S 


2002 Aug 4 


2685 


13.7 


Chandra/ACIS-S 


2002 Sep 9 


2683 


14.3 


Chandra/HRC-S 


2002 Nov 8 


2797 


49.4 


Chandra/HRC-S 


2006 May 27 


6769 


41.1 


Chandra/ACIS-S 


2008 Aug 7 


9132 


144.1 


Chandra/ACIS-S 


2008 Aug 10 


9133 


55.2 



Due to the relatively short exposure, just six sources had 
sufficient counts for a spectral fitting analysis, which in- 
clu ded the original globu lar cluster MSP, PSR B1821- 
24. iRutledge et al.l (|2004[ ) have used the superb timing 
capabilities of the HRC-S instrument on board Chan- 
dra to study the X-ray pulsations from B1821-24. More 
recently, this MSP was again targeted with HRC-S to 
calibrate the absolute timing precision of this detector 
(unpublished data). 

In this paper, we present new deep Chandra X-ray Ob- 
servatory ACIS-S observations of M28 and its 12 known 
pulsars. This investigation sheds further light onto the 
general X-ray properties of the Galactic population of 
MSPs. A comprehensive analysis of the other X-ray 
source classes in M28 will be presented elsewhere (M. van 
den Berg, et al., in preparation). The present work is out- 
lined as follows. In §2, we describe the observations and 
data reduction and analysis procedures. In §3 we present 
the imaging analysis and detection methods of the X-ray 
counterparts to the MSPs. In §4, we conduct a spec- 
troscopic and timing study of PSR B1821-24 (M28 A), 
while in §5 we focus on the peculiar binary PSR J1824- 
2452II (M28 H). In §6 we discuss the remaining pulsars 
in the cluster. We offer conclusions in §7. 

2. OBSERVATION AND DATA REDUCTION 

The new Chandra dataset presented herein was ac- 
quired during two separate observations in Cycle 9, on 
2008 August 7 (ObsID 9132) and 2008 August 10 (Ob- 
sID 9133) for 144.1 and 55.2 ks, respectively. For both 
exposures, the optical center of the cluster at a = 
18^24 '"32.89" and S = -24°52'11.4" (jShawl fc White! 
[19861) was positioned 1' off-axis from the nominal aim 
point of the back-illuminated ACIS-S3 CCD, configured 
in VFAINT telemetry mode. We also include the three 
Cycle 3 ACIS-S observations of 14.3, 12.9, and 13.7 
ks (ObsIDs 2683, 2684, and 2685, respectively) from 
IBecker et al] ()200^. also taken in VFAINT mode. In the 
analysis of PSR B1821-24 (see §4) we make use of the 
two HRC-S exposures of 49.4 and 41.1 ks (ObsIDs 2797 
and 6769, respectively) in the special SI mode, which 
permits a timing resolution of ~16 /is. All observations 
used herein are summarized in Table 1. 

The data re-processing, reduction, and analysis were 
performed using CIACH 4.1.2. Starting from the 
level 1 data products of the ACIS-S observations, we 
first removed pixel randomization from the standard 
pipeline processing to aid in source disentanglement in 
the crowded cluster core. In addition, for the pur- 

C handra Interactive Analysis of Observations, available at 
|http : //c xc .harvard. edu/ciao/ 



poses of faint source detection we applied the background 
cleaning algorithm specific to the VFAINT telemetry 
mode. However, this procedure tends to reject real source 
counts for relatively bright sources. Consequently, the 
background-cleaned data were not used for the spectral 
analyses discussed in §4 and 5. Before coadding the 
five ACIS-S images, the individual aspect solutions were 
reprojected relative to the longest observation (ObsID 
9132) using the brightest sources in the cluster to cor- 
rect for any differences in the absolute astrometry be- 
tween the individual observations. The summed ACIS-S 
data results in a 237.1 ks effective exposure time, one of 
the deepest X-ray observations of a globular cluster, with 
a limiting sensitivity of Lx '~ 1 x 10^° ergs s~^ for a 1 
keV source at the distance of I? = 5.5 kpc and hydrogen 
column density iVn = 2.4 x 10^^ cm~^. 

For the spectroscopy and variability analyses of the X- 
ray counterparts to the MSPs, we extracted the emission 
from polygonal (roughly circular) regions that enclose 
'--^90% of the total source energy at 1.5 keV. To permit 
spectral fitting in XSPECFl. the extracted source counts 
in the 0.3-8 keV range were grouped in energy bins so 
as to ensure at least 15 counts per bin. The background 
was taken from three source-free regions in the image 
around the cluster core. For the pulse profile of pulsar A 
and variability analysis for pulsar H, the photon arrival 
times were reduced to the solar system barycenter, using 
the CIAO tool axbary assuming the DE405 JPL solar 
system ephemeris. 

3. IMAGING ANALYSIS AND SOURCE DETECTION 

Figure 1 shows the co- added image of all ACIS-S ex- 
posures (237.1 ks in total) of the core of M28 in the 
0.5-2 keV band (where most of the emission from MSPs 
is expected), with 1" circles centered on the 12 current 
best radio pulsar timing positions (S. Begin, et al. in 
preparation). It is apparent that the unresolved diffuse 
em ission in the s hallow er ACIS-S observations reported 
by IBecker et al.l ()2003[ ) is mostly resolved into multiple 
faint sources, with some apparently coinciding with the 
known radio pulsars in M28. To formally identify the X- 
ray counterparts to the pulsars, we used the CIAO tool 
wavdetect, which correlates the image with Mexican- 
hat wavelets of various scales to identify sources. We 
also used the PWDetectP^ script (Damiani et al. 19973BI), 
which is generally more effective at identifying faint 
sources near much brighter sources. Subs equently, we 
empl oyed the IDL script acis.extraclQ (jBroos et al.l 
l2010f ) to confirm the validity of the source detections re- 
ported by wavdetect and PWDetect and refine the source 
positions. 

For PSR B1821~24 (henceforth MSP A, M28 A, or pul- 
sar A) we find an offse t between the radio pulsar position 
from lRav et al.l (|2008[ ) and the known X-ray counterpart 
of Aa = —0.15" and AS = —0.26" in right ascension and 
declination, respectively. Thus we apply this correction 
to the astrometric frame of the X-ray dataset, which re- 
sults in significantly better matches between most of the 
radio MSPs and X-ray sources. The positional uncer- 

Available at http : //heasar c . nasa . gov/docs/ xanadu/xspec/index ■ html | 
^■^ Sec http : //www . astropa . unipa . it/progetti_ricerca/PWDetect/ 
for details. 

Available at |http : //www ■ astro . psu . edu/xray/docs/TARA/ae_user s.guideT 
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Fig. 1. — (Left) Chandra ACIS-S coadded 0.5—2 keV image of the globular cluster M28 with the positions of the 12 known radio pulsars 
in the cluster labeled. (Right) The inner core of M28, with 1" radius red circles centered on the radio pulsar timing positions and blue 
polygons representing the extraction regions of the nearest X-ray source to each pulsar, with the exception of pulsar G for which the nearest 
three sources are shown. The magenta ellipses show the TEMPO radio timing position uncertaint for pulsars D and I, multiplied by a 
factor of 2 and 5, respectively (see text for details). For all other pulsars, the radio position uncertainties are much smaller (<0.1"). The 
image greyscalc corresponds to number of counts increasing logarithmically from (white) to 2553 (black). North is up and east is to the 
left. 



tainties of the detected X-ray s ources were compu ted us- 
ing the empirical relation from lHong et al.l ()2005[ ). which 
yields the 95% error radius, ^95%. We note that the un- 
certainties in the radio timing positions for the M28 pul- 
sars are generally ^0.1". Therefore, the dominant source 
of uncertainty in the X-ray counterpart matches is the 
error in the X-ray source position, with the exception of 
pulsars D and I. Pulsar D is significantly younger (based 

on its characteristic age r = P/2P « 26 Myr) then the 
other pulsars in M28 and has a comparatively long spin 
period of P = 79 ms. As such, it exhibits more timing 
noise, resulting in an uncertainty in its radio position, es- 
pecially in declination, comparable to rg5% (^0.5" — !")• 
The error ellipse shown in Figure 1 was obtained by dou- 
bling the formal uncertainties in the pulsar position ob- 
tained with TEMPO. In the case of MSP I, the presence 
of random and irregular radio eclipses, apparent long- 
term variations in the orbital motion, and infrequent de- 
tections of the pulsar pose major difficulties in deriv- 
ing a unique phase-connected timing solution and thus 
a precise radio timing position (see S. Begin et al. in 
preparation, for details). Thus, the error ellipse shown 
in Figure 1, taken to be five times the formal position 
uncertainty derived from a TEMPO fit with a xt over 
7, represents the current best estimate of the location of 
pulsar I, which may be subject to significant change. We 
note that the best timing fit of the eclipsing binary MSP 
H also has a large xtj implying that the uncertainties in 
the position reported by TEMPO are not realistic esti- 
mates of the true errors. Fortunately, the recent identifi- 
cation of its optical counterpart by Pallanca et al. (2010) 
provides a much more reliable localization of M28H. 

As evident from Figure 1 and Table 2, MSPs A, C, 
H, J, and K are found well within the 95% error cir- 



cles of X-ray sources, with differences between the X-ray 
and radio positions of <0.24". Thus we deem the X-ray 
identifications of these pulsars as secure. For pulsar D, 
a single X-ray source falls within the radio timing error 
ellipse. As the source is found away from the crowded 
cluster core, it is most likely the X-ray counterpart to 
the pulsar. Pulsar F is located just 0.04" outside of rg5% 
of the sole X-ray source in the vicinity, suggesting that 
the two are related. It is likely that the small positional 
discrepancy is due to a slight rotational offset between 
the X-ray and radio astrometric frames. However, as 
pulsar A is the only bright X-ray source with a radio 
counterpart in this field, it is not possible to unambigu- 
ously correct for this effect. Pulsar E is found 0.12" out- 
side of the 95% error circle of the nearest X-ray source. 
As this MSP is located in the crowded cluster core, this 
association is less certain. Although faint emission is 
coincident with the position of MSP B, the source detec- 
tion algorithms fail to identify an X-ray counterpart to 
this pulsar. The detection of MSP I is rendered difficult 
by its position between MSP K and one of the b right- 
est sources in the cluster (23 in lBecker et al.ll2003[ ). The 
centroid of the X-ray counterpart to pulsar K is nearly 
coincident with its radio position (|Pr — Px\ =0.03"), 
implying that pulsar I is either significantly fainter than 
K or is not at the location derived from radio timing. 
MSP L is nearly coinci dent with a mo derately luminous 
X-ray source (22 iu iBecker et al.|[2003l) . As discussed in 
§6, however, based on the properties of this pulsar as in- 
ferred from radio observations and the temporal behavior 
of the X-ray emission from source 22, the two objects are 
most likely not physically associated. For MSP G, three 
X-ray sources are found ~1" away (see Figure 1), mak- 
ing the X-ray-radio association ambiguous as the X-ray 



4 Bogdanov et al. 



TABLE 2 

Net Counts and Luminosities for the pulsars in M28 



PSR J1824-2452 


Soft'' 
(0.3-1.5 keV) 


Hard'' 
(1.5-6 keV) 


(0.5-8 keV) 


'•95%'' 
(") 


\Pr - Pxr 
(") 


Detection 


A 


2106 


3800 


1147.0 


0.28 


0.0 


Yes 


B 






<2 






No 


C 


11 


6 


1.7 


0.39 


0.21 


Yes 


D 


8 


2 


1.0 


0.44 


0.99 


Yes 


E 


8 


8 


2.0 


0.40 


0.49 


7 


F 


13 


10 


1.2 


0.50 


0.54 


Yes 


Qf 


27 


6 


<10.8 


0.33 


1.19? 


7 


H 


23 


49 


16.9 


0.33 


0.10 


Yes 










$1 






No 


J 


37 


21 


1.3 


0.34 


0.20 


Yes 


K 


26 


20 


5.4 


0.35 


0.03 


Yes 




350 


536 


<189.2 


0.30 


0.43 


No 



^ Background-subtracted net counts extracted from regions enclosing ~90% of the total energy at 1.5 keV from the nearest X-ray source 
to each pulsar. For MSPs I and B, for which no X-ray counterpart is detected, the source counts were extracted from a circle of radius 
1.7" centered on the radio timing position. 

^ Estimated unabsorbed X-ray luminosity in units of 10^" ergs s~^ in the 0.5-8 keV band based on Cash statistic (Cash 1979) fits to the 
spectra with a power-law model for Nn = 2.4 X 10^^ cm~^, and D = 5.5 kpc. For MSPs A, G, H, and L the values were derived from 
spectral fitting using the statistic with A^h as a free parameter (see Table 3). 

Positional uncertainty (95% confidence error circle radius) in arcseconds of the putative X-ray counterpart to each radio MSP, determined 
using Equation 5 from Hong et al. (2005). For pulsars B and I there are no plausible X-ray counterparts in the vicinity. 

Angular separation in arcseconds between the radio timing position of each MSP (from S. Begin, et al., in preparation) and the centroid 
of the nearest X-ray source. For pulsar A, the difference is zero as this pulsar was used to align the radio and X-ray astrometric frames. 
° Emission strongly contaminated by neighboring bright source. 
' Possibly an unresolved blend of two or more unrelated sources. 



photons coincident with its radio position may be due to 
collective emission from the three nearby X-ray sources. 
It is also difficult to reliably constrain the intrinsic pul- 
sar flux as there is no appropriate region from which to 
extract a representative background. Consequently, in 
order to place a crude upper limit on the flux, we ex- 
tracted source counts within 1" of the radio position. 

Given the high X-ray source number density in the 
cluster core, it is important to investigate the possibility 
that X-ray sources coincide with most of the radio pulsar 
positions purely by chance (with the exception of M28 A, 
for which the radio-X-ray association is well established). 
We have applied systematic offsets in the range 0.5— 5" in 
all directions of the radio pulsar positions with respect to 
the X-ray dataset. We find that any offset consistently 
results in two or less source matches. The same holds 
true for all angles of rotation of the radio pulsar positions 
about the cluster center or about the position of pulsar A. 
This finding suggests that the majority of the radio-X- 
ray source matches are indeed genuine and is consistent 
with pulsars L and/or G being possible false matches. 

Table 2 summarizes the background-subtracted X-ray 
net source counts extracted from polygonal regions that 
enclose ^^90% of the total energy at 1.5 keV for the pu- 
tative X-ray counterparts of the pulsars. For MSPs I 
and B, for which no X-ray counterparts are detected, we 
extracted source counts from circles of radius 1.7" (en- 
closing 95% of the energy at 1.5 keV) centered on the 
radio timing positions. The background event contribu- 
tion was determined from 1.7" circles that are th e same 
radial distance from source 23 ()Becker et al.ll2003D as the 
MSPs. In both cases the net source counts are consistent 
with zero. 

For the X-ray-faint pulsars (C, D, E, F, J, and K), 
we estimate the unabsorbed X-ray luminosity by fitting 
a powerlaw spectrum with Nh = 2.4 x 10^^ cm^^ us- 



ing the Cash statistic (jCashI 119791 ) in XSPEC. This par- 
ticular spectral model was chosen because it reproduces 
the rough spectral shape of a two-temperatu re thermal 
model observed in nearby MSPs (Becker & A schenbachI 
|200iLZavlin 2006; Bogd anov fc Grindl av 2009) quite weh 
for crude photon statistics, with spectral photon index in 
the range 2 — 3.5. In the case of MSPs A, G, H, and L 
the luminosities were derived from the spectral fits dis- 
cussed in §4 and 5, although for G and L these values 
most probably do not correspond to the intrinsic pulsar 
values as a consequence of source confusion/blending. 

We note that in the combined Chandra HRC-S data 
(totaling 90 ks), of the 12 pulsars only pulsar A is de- 
tected, owing to the ~3 times lower sensitivity at 0.5-3 
keV and the higher background of HRC compared to 
ACIS. 

4. PSR B1821-24 

PSR B1821-2 4 is the first MSP to be discovered in a 
globular cluster ()Lvne et al.lll987f ). This solitary pulsar 

with P = 3.05 ms, P = dP/dt = 1.61 x lO'^^ s s-\ is 
~100 times younger than most "recycled" MSPs, based 
on its characteristic age (tc — P/2P = 30 Myr). It is 
also the most energetic MSP known, with spin-down lu- 
minosity E = 2.2 X 10^^ ergs s^^ and X-ray luminosity 
Lx ~ 1.3 X lO^'^ ergs s ~^ (phase integrated for 0.5-8 keV, 
see lBecker et al.ll2003D . 2-3 orders of magnitude greater 
than typical MSPs. Only two other X-ray detected 
MSPs, PSRs B1937+21 and J0218-H 4232 (both field 
MSPs) , have comparable values of Lx (jCusumano et al.] 
[2003; Webb et al. 2004). 

Because of the bright nature of M28 A, --4% of the X- 
ray photons from this source in the ACIS-S observations 
are expected to suffer from pileup (based on PIMM^3)- 

See |http : / /cxc . harvard . edu/toolklt/plnms ■ j sp . | 



X-rays from Millisecond Pulsars in M28 



5 



0.01 



10-- 



PSRB1821-24 




^^iait,,tiii,,i,.j,Miy.JjyiiiiyiiM]kk 



0.5 



1 2 
Photon energy (keV) 



Fig. 2. — Phase-integrated Chandra ACIS-S X-ray spectrum of 
PSR B1821— 24 fitted with an absorbed power-law spectrum. The 
lower panel shows the best fit residuals. See text and Table 3 for 
best fit parameters. 

Pileup occurs when two or more photons arrive within 
the same detector pixel within a single ACIS frame in- 
tegration time and are registered as a single event. The 
detected energy of such an event is approximately equal 
to the sum of the energies of the individual photon events 
of which it is comprised. The result is a distortion of the 
intrinsic source spectrum (Davis 2001). To ascertain the 
severity of this effect we have included the pileup model 
provided in XSPEC in our spectral fits. We find that 
pileup has negligible impact on the spectrum of PSR 
B1821-24, as the best fit values of the model parame- 
ters with and without pileup are well within la of one 
another. Based on this, in what follows, we ignore the 
effect of pileup. 

As expected, the phase-integrated X-ray spectrum of 
PSR B1821-24 is very well described by a pure absorbed 
power- law (Figure 2). The best fit parameters are F = 



(2.55i;j i^) X 10 cm-2 with x 



1.226t^:°^^, 

0.93 for 253 degrees of freedom. The uncertainties quoted 
are la for one interesting parameter. The unabsorbed 
energy flux is Fx = (S.SOt^ J^) x lO'^^ g^gg g-i (q 3, 
8 keV), corresponding to a luminosity of Lx = 1.4 x 
10'^^&{D/5.5 kpc)^ ergs s~^, where & is the (unknown) 
portion of the sky covered by the X-ray emission beams. 

By analogy with most MSPs, M28 A may have heated 
magnetic polar caps. As a result, low-level thermal emis- 
sion could be present in the predominantly non-thermal 
spectrum of PSR B1821-24. Introducing either a black- 
body or hydrogen atm osphere model (NSATMOS from 
iMcClintock et all 120041 ) with temperatures in the range 
(0.5 — 3) X 10^ K yields acceptable fits, although given 
the high-quality of the pure power-law fit, the addition 
of the thermal component is not warranted by the data 
(as confirmed by an F-test). This allows us to set a firm 
upper bound on any thermal component of ^-^1% of the 
total energy flux in the 0.3-8 keV range. The implied 
thermal luminosity of 540'^^ ergs s~^ is comparable to 
most nearby MSPs (|Zavlinl 12006: Bogdanov fc Gr indlavl 
[200l and those in 47 Tuc (Bogdanov et al. 20061). On 
the other hand, in some theoret ical pulsar models (see, 
e.g., iHarding fc Muslimovl 120021 ) . the energetics of the 
magnetosphere of PSR B1821-24 may not be conducive 
to the production of a substantial return current of par- 
ticles to heat the polar cap surface so this pulsar may in 



f act have much cooler polar caps. 

iBecker et al.l ()2003D have reported intriguing but 
marginal evidence (at 98% confldence) for a broad emis- 
sion feature centered at 3.3 keV, which they interpreted 
as cyclotron emission from the magnetosphere above the 
pulsar polar cap if the magnetic field strongly deviates 
from a centered dipole. We have reprocessed the archival 
ACIS-S data using the latest Chandra X-ray Center cal- 
ibration products and find no evidence of this feature. 
As seen from Figure 2, the broad spectral feature is ab- 
sent in the total ACIS-S spectrum as well. Therefore, 
we attribute the deviation from a pure power-la.w of the 
PSR B1821-24 spectrum found bv lBecker etldl (|2003D to 
an unmodelled instrumental feature that has since been 
incorporated into the telescope calibration model. 

4.1. The Unpulsed X-ray Emission from PSR B1821-24 

The bulk of X-rays from PSR B18 21-24 (^85%) 
are e mitted in two very narrow pulses (jRutledge et al.l 
|2004| ). indicative of highly beamed, non-thermal ra- 
diation from the pulsar magnetosphere, in stark 
contrast with t he broad, s o ft thermal pulsations 
of most MSPs (iZavlinI [200l iBogdanov et al l l200l 



iCameron et"alll2007t iBogdanov fc Grindlavl 120091 ) . Sur- 
prisinglv. iRutledge et al.l (|2004l ) have found an appre- 
ciable unpulsed component, contributing ~15% of total 
photon flux. 

Using the two HRC-S timing observations of PSR 
B 182 1-24 we can place constraints on the origin of the 
puzzling unpulsed emission. To this end, we combined 
the pulse proflles from the two HRC datasets as fol- 
lows. We flrst folded the barycentered arrival times of 
the events extracted from a 1" region centered on the 
radio pulsar in TEMPC(3 using the timing ephemeris 
from iRav et al.l ()2008i) for each dataset separately. Due 
to the ~3. 5-year gap between the two HRC-S observa- 
tions it is not possible to phase-connect the two datasets 
accurately using the available radio timing ephemerides. 
Thus, to determine the relative phase misalignment, we 
cross-correlated the folded pulse profiles from both ob- 
servations. The individual and combined pulse profiles 
are shown in Figure 3. The summed HRC-S pulse pro- 
file yields 706 photons, a 1.8- fold improvemen t in ph oton 
statistics over the data from IRutledge et al.l (|2004[ ). Of 
these, 16 ± 0.4 photo ns are due to background. As in 
IRutledge et~all I12004D . the unpulsed (DC) level in the 
summed profile was determined by considering the emis- 
sion ±0.15 in phase away from the peaks of the two pulsas 
to minimize the contribution from the pulsed component. 
After subtracting the background contribution, we find 
that 17.5% ± 4% of the source photons are unpulsed. 
Varying the choice of phase intervals used to determine 
the unpulsed level results in a spread of ±1% around the 
nominal value. 

Although the HRC has poor intrinsic energy resolu- 
tion, in principle, it is still possible to determine any 
gross difference between the spectra of the pulsed and 
unpulsed emission using hardness ratioJ^. The hard- 
ness ratios, however, do not show any substantial differ- 
ences in the spectra of the pulsed and unpulsed portions. 
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Fig. 3.— Chandra HRC-S pulse profiles of PSR B1821-24 from 
ObsIDs 2797 and 6769 (top panel) and the summed pulse pro- 
file (middle panel). The bottom panel shows the summed profile 
smoothed using a 4-point moving average. Note the asymmetry of 
the wings of the two pulses. The choice of phase zero is arbitrary. 
Two cycles are shown for clarity. 



implying that the unpulsed emission is not drastically 
softer or harder than the pulsed emission. Nonethe- 
less, the high-quality phase-integrated ACIS-S contin- 
uum provides stringent constraints on the spectral prop- 
erties of the unpulsed emission, in particular, if it is of 
a thermal nature or not. The best-fit absorbed power- 
law model derived from the ACIS-S spectrum, convolved 
with the HRC-S effective area curve, reveals that the bulk 
of photons from PSR B1821-24 in the HRC-S datascts 
are from the 0.4-2 keV range. We fitted the ACIS-S 
spectrum with additional thermal components to deter- 
mine their photon flux contribution in this energy range. 
Converting the resulting values to the equivalent HRC-S 
count rates indicates that any thermal emission can con- 
tribute no more than ^6% of the total number of photons 
from M28 A in the HRC-S observations. This leaves at 
least ~7-15% of the photons unaccounted for. 

It is interesting to speculate on the origin of this excess 
unpulsed emission. In many young, energetic pulsars, 
with the Cra b pulsar being the most notable example 
(|Kirsch et al.|[200Q) . the high -energy non-thermal pulse 
profiles often exhibit a "bridge" of emission connecting 
two pulses, suggesting that the two pulses originate from 
the same emission cone. However, for PSR B1821-24 
the count rates in the regions between the two pulses are 
consistent with being identical. The apparent lack of ex- 
cess of photons in one of the regions implies that bridge 
emission is not responsible for the seemingly unpulsed 
photons. It should be noted however that both pulses 
exhibit tails of emission from their trailing edges (see 
bottom panel of Figure 4) , which could plausibly extend 
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^ Spectrum likely dominated by unrelated X-ray source(s). 
^ Hydrogen column density along the line of sight in units of 
10^^ cm-^ 

" Unabsorbed flux in units of 10 ergs cm ^ s ^ in the 0.3-8 
keV band. 

to cover most of the rotational period. Indeed, certain 
high-energy emission models of pulsars, such as the two- 
pole caustic model, predict low-level non-thermal magne- 
tospheric emission spanning the entire period of rotation 
(jVenter et al.|[2010f ). 

A plausible alternative origin for the unpulsed photons 
is synchrotron emission from an unresolved wind neb- 
ula surrounding this rather energetic pulsar. By anal- 
ogy with pulsars with c omparable spin-down luminosi- 
ties (see iKargaltsev fc Pavlov..2008t , for a review), the ex- 
pected angular size of a wind nebula surrounding pulsar 
M28 A, assuming a distance of 5.5 kpc, would be 
Such a nebula would be difficult to resolve, especially 
considering the source crowding in the cluster core. If 
the unpulsed emission arises solely due to the wind neb- 
ula, its implied X-ray luminosity is ipwN ~ 10'^^ ergs 
s~^, similar to what is observed for other energetically 
comparable pulsars ( see in particular Tables 1 and 2 in 
iKargaltsev fc Pavlov|[2008l ). 

Conceivably, the unpulsed emission could be caused 
by small-angle scattering of X-rays p assing through dus t 
grains in the interstellar medium (jSmith et al.l 1200^ . 
The additional path length of the scattered X-rays would 
induce a time delay in the photon time of arrival, caus- 
ing them to become out of phase relative to the unscat- 
tered photons and thus form an unpulsed component. 
This effect could also account for the trailing tails of 
the two pulses. The location of M28 in the Galactic 
disk {I — 7.8°, b = 5.6°) suggests that dust scatter- 
ing is not neglegible. We estimate the optical depth 
to scattering using the empirical relati on Tisra.t.ikfiV = 
At/(0 .079 ± 0.003) - (0.052 ± 0.019) (.Predehl fc SchmittJ 
Il995f ). where Ay is the visual interstellar ext inction. For 
Ay ~ 1.4 mag towards M28 ()Davidgelll996l) . this yields 
T r-^ 0.05 implying that I/Io = e""^ ~ 5% of the X-ray 
photons undergo scattering. 

Finally, it is possible that at least a portion of the 
unpulsed component from PSR B1821-24 is due to one 
or more fainter, unrelated X-ray sources within ~1" of 
PSR B1821-24. However, using wavdetect we find that 
the centroid of the unpulsed emission is offset by only 
0.07" from the centroid of the total emission, suggesting 
that source contamination is unlikely. 



5. PSR J1824-2452H 

Apart from PSR B1821-24, the binary PSR J1824- 
2452H is the only other unambiguously detected MSP 
with sufficient counts for a spectral fitting analysis. PSR 
J1824-2458H is in a 10.4-hour circular orbit around 



X-rays from Millisecond Pulsars in M28 



7 



2x10-' 



5x10-5 



2x10-5 
2 

X 
-2 



2x10-' 
^ 10-' 
U 5x10-* 



PSR 1824-2452G 



PSR 1824-2452H 



PSR 1824-2452L 




^/ llt| t# i f ll % i ++ 



0.5 



1 2 
Photon energy (keV) 



Fig. 4.— X-ray spectra of PSR J1824-2452H and of the X-ray 
sources nearest to PSRs J1824-2452G and L, fitted with power-law 
spectra. The bottom panels show the best fit residuals. See Table 
3 for best fit parameters. 



what appears to be a non-degenerate star (jPallanca et al.l 
l2010f ) and exhibits highly variable and irregular radio 
eclipses (S. Begin et al. in preparation). This atyp- 
ical MSP binary system is thus likely the result of a 
past dynamic al binary exchange enc ounter in the dense 
cluster core (jCamilo & Rasid 120051), like PSR J1740- 
5340 in NGC 6397 ( D'Amico et al. 2001; Grindlay^al] 
I2002t iBogdan ov eFal.ii20ia IHuang fc Becker 2010) and 
PSR J0024-7204W i n 47 Tuc (ICamilo e t al. 2001 
lEdmonds etHI 120021: IBogdanov et aTTMosI) . Alterna- 
tively, it may have just emerged from the low-mass 
X-ray binary phase with the MSP re-activating as a 
rotation-powered pulsar, as is also possibly the case for 
PSR J1023-h0038 (see Archibald et al. 2009). The X- 
ray properties of PSR J1824-2452H are thus of particu- 
lar int erest as these systems show unusual X-ray prop- 
erties (IStappers et alj|2003l: IBogdanov et al.|[2005ll2010l: 
I Archibald et al.ll2010t) . We fitted the spectrum of this 
MSP with a power-law model (middle panel of Figure 4), 
yielding TVh = (2.ltl:i) x lO^^ cm'^, F = imtoil, and 
unabsorbed flux of Fx = i'i-Stl^) x 10'^^ (0.3-8 keV), 
with xl = 0.87 for 5 de grees of freedom. The nominal 



Time (hours) 
10 



20 



0.0008 



o 0.0006 



0.0004 



c 

3 

o 
u 



0.0002 



PSR J1824-2452H 



0.3-8 keV 



_l I I I I I L_ 



0.5 1 1.5 

Orbital phase 



Fig. 5.— X-ray lightcurve of PSR J1824-2452H folded at the 
iO. 4-hour binary period of the pulsar. The orbital phase is defined 
based on the radio pulsar timing convention in which superior con- 
junction occurs at (/) = 0.25. Two orbital cycles are shown for 
clarity. 

best fit flux translates to a luminosity of Lx = 1-7 x 10^""^ 
ergs s~^ for a 5.5 kpc distance. Despite the limited 
photon statistics the spectrum is poorly reproduced by 
a single-temperature thermal model. Applying models 
with two or more components offers no meaningful con- 
straints on any model parameters. 

For this binary MSP we can also investigate any large 
amplitude X-ray variability as a function of orbital phase, 
which could potentially provide clues regarding the ori- 
gin of the X-ray emission. We have used the radio tim- 
ing ephemeris of PSR J1824-2452H (S. Begin et al., in 
preparation) to fold the X-ray dataset at the binary pe- 
riod. The resulting lightcurve (Figure 5) shows com- 
pelling evidence for variability. For four and five phase 
bins over the binary orbit, the variability is significant at 
a 4cr and 2.7a level, respectively, based on a test. A 
Kuiper test ([Paltani 2004) on the unbinned lightcurve, 
weighted to account for the non-uniform exposure across 
the orbit, indicates a 5.8 x 10^* {3.2a) probability that 
photons being drawn from a constant distribution would 
exhibit this level of non-uniformity. It is interesting 
to note that the minimum around (f) — 0.25 in the X- 
ray lightcurve is roughly coincident with the position 
of the radio eclipse. Moreover, the rough shape of the 
lightcurve is remarkably similar to those of the analogou s 
binaries PS Rs J0024-7204W (IBogdanov et all [2005h . 
J1740-5340 (IBogdanov et al.l |2010|), and J1023-h0038 
(|Archibald et al.r2010[ ). Thus, the peculiar nature of 
the PSR J1824-2452H binary system, the implied X-ray 
luminosity, relatively hard spectrum compared to most 
MSPs observed in X-rays, and the suggested variabil- 
ity at the binary period point to an intrabinary shock 
origin of the X-rays associated with this object. The 
minimum in the X-ray light curve at (/) « 0.25 can be 
plausibly interpreted as a geometric occultation (either 
partial or total) of the shock by the secondary star. Al- 
ternatively, it could arise due to relativistic beaming of 
the X-rays from the radiating particles in the shock. Due 
to the limited photon statistics and the number of or- 
bits covered ('^6) it is difficult to establish whether the 
fiux modulation is truly periodic and stable over many 
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orbit s As seen in PSR J0 024-7204W (|Bogdanov et alj 
l2005t iCameron et al.|[2"007[ ). the orbital variation of the 
X-ray flux may actually evolve on timescales spanning 
from months to years. 

6. OTHER MSPS 

The ^10 times larger absorbing column towards M28 
compared to 47 Tuc and '^1.2 times greater distance 
make it more diffic ult to study soft X-r ay sources such 
as thermal MSPs (Bog danov et al.l 12006*) in a compara- 
ble exposure. As a result, most of the detected pulsars in 
M28 do not have adequate photon statistics for detailed 
spectral analyses. Nonetheless, we can still constrain the 
origin of the observed X-rays by considering the summed 
spectrum of pulsars B, C, D, E, F, J, and K (shown 
in Figure 6). Indeed, the combined X-ray continuum 
of these seven pulsars is poorly described by a single- 
temperature thermal spectrum, while a pure power-law 
model requires A'h > 4 x 10^^ cm for an acceptable fit, 
substantially higher than the value suggested by the fit to 
the spectrum of PSR B1821-24. A composite blackbody 
plus power-law model for N — 2.4 x 10^^ cm~^ gives an 
acceptable fit with kT = 0.26±0.03 keV, i?cff = 0.19t!]:ij^ 
km, F = 0.8 and xl = 0.86 for 10 degrees of free- 
dom. Note that a power-law component does not nec- 
essarily imply non-thermal emission since the superposi- 
tion of several multi-temperature thermal MSP spectra 
can produce this s pectral shape for lim ited photon statis- 
tics (|Zavhnll2006t iBogdanov fc Grind lav 2009). The de- 
rived unabsorbed flux of 7.5 x 10~^^ ergs cni^^ s~^ (0.3-8 
keV) implies a mean X-ray luminosity of 3.9 x 10^° ergs 
s^^ for the seven pulsars or 4.5 x lO'^" ergs s~^ if we ex- 
clude pulsar B. For comparison, the mean unabsorbed 
X-ray luminosity of the 11 thermal MSPs in 47 Tuc with 
unconfus ed and uncontan iinated spectra is ^ 3.2 x lO'^'^ 
ergs s~^ (jBogdanov et al.l[200QV This suggests that the 
MSP populations of the two clusters are quite similar in 
terms of X-ray properties. 

The non-detection of pulsars B and I implies that 
these pulsars may be below the sensitivity of the total 
ACIS-S exposure, with Lx < I x 10^° ergs s'^ (0.3- 
8 keV), perhaps similar to the unusually sub-luminous 
PSR B1257-hl2 wit h Lx ^ 2.5 x 10^ ^ ergs s^^ in the 
same energy band UPavlov et al.ll2007D . The low lumi- 
nosity of the non-recycled pulsar D is simi lar to those 
of co mparably old nearby pulsars (see ,Kar galtsev et al.l 
l2006[) . 

For pulsars G and L, there are sufficient counts within 
1" to conduct spectral fitting and search for variability, 
although as noted in §3, the association of this emis- 
sion with the pulsars is questionable. For both sources, 
the spectra are well described by an absorbed power-law 
model (see Table 3 and Figure 4). In the case of pulsar G, 
there is no statistically significant variability. As for PSR 
J1824-2452L, the X-rays extracted from within 1" of the 
radio timing position exhibit variability on time-scales of 
order hours. In addition, there is a substantial increase 
(a factor of ^3) in the count rate from the three 2002 
ACIS-S observations to the two observations in 2008. In 
principle, such large-amplitude variability from an MSP 
binary could arise due to interaction of the pulsar wind 
with gas from the companion star. However, this MSP 
is not unusual and exhibits no radio eclipses, which are 
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Fig. 6.— Summed X-ray spectrum of PSRs J1824-2452B, C, D, 
E, F, J, and K fitted with a blackbody plus powerlaw model. The 
bottom panels show the best fit residuals. See text for best fit 
parameters. 

associated with intrabinary shocks in all known cases. 
Therefore, the variable X-rays are most probably from 
source 22 from lBecker et ahl i|2003l ) . which is not related 
to the pulsar. 

7. CONCLUSION 

We have presented deep Chandra observations of the 
sample of known pulsars in the globular cluster M28. 
Seven of the 12 known pulsars in the cluster are se- 
curely detected (A, C, D, F, H, J, and K) two are possi- 
bly detected (E and G), while three are not detected 
(B, L, and I). Most of the MSPs firmly or possibly 
detected appear to be relatively faint and soft sources 
most likely due to thermal polar cap emission. In this 
sense, they appear to be akin to the majority of "recy- 
cled" MSPs in 47 Tuc (IBogdanov et al.l'l2006l) and the 
field of the Galaxy ()Zavlinll2006l: IBogdanov fc Grindlav! 
120091 ) ■ with soft spectra and Lx ~ 4 x lO^o-si ergs 
s~-^ (0.3-8 keV). This provides further evidence that 
in typical MSPs, both in globular clusters and the 
field of the Galaxy, the X-ray emiss ion is primarily 
of a t hermal nature. As sh o wn bv | Pavlov fc ZavlinI 
(IT9971. Za ^lin fc Pa vlov' [TME). Bog danov et al.l (I2007D . 
and Bogd anov et al . (2008), realistic modelling of the 
X-ray spectra and pulsations of these objects may po- 
tentially allow stringent constraints on the neutron star 
equation of state through measurement of the mass-to- 
radius ratio of these stars. As such, MSPs warrant further 
investigation in the X-ray band. 

The most detailed X-ray spectrum to date of the en- 
ergetic PSR B1821-24 is found to be well described by a 
pure power-law with F — 1.2, with no requirements for 
additional components. We set a limit on any thermal 
X-ray emission from the MSP of <10'^^ ergs s~ for 0.3-8 
keV. The greatly improved ACIS-S spectrum of this pul- 
sar shows no evidence for the previously reported spectral 
feature at ^3.3 kcV. A re-examination of the X-ray pulse 
profile of PSR B1821-24 reveals that the observed un- 
pulsed emission is unlikely to be of thermal origin and is 
due to either low-level non-thermal magnetospheric emis- 
sion, unresolved nebular emission, and/or dust scattering 
of the pulsed X-rays. 

PSR J1824-2452H appears to share common X-ray 
characteristics (namely, variable and hard non-thermal 
emission with Lx ~ 10^^ ergs s~^) with PSRs J0024- 
7204W, J1023-f0039, and J1740-5340. Presumably, 
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in these peculiar binary systems the X-rays originate 
from an intrabinary (synchrotron) X-ray emitting shock 
formed by interaction of the MSP wind and matter 
from the secondary star (jArons fc TavanilfTool . Fur- 
thermore, this class of objects exhibits remarkably sim- 
ilar X-ray and optic al properties to the accreting X-ray 
MSPs in qu iescence (jCampana et al.ll2004l : iHeinke et alj 
120071 12009() and thus can potentially offer vital clues into 
the little understood activation mechanism for rotation- 
powered MSPs. 



We would like to thank Z. Medin, A. Harding, and 
A. Archibald for useful discussions and the referee 
F. Camilo for comments. The work presented was funded 
in part by NASA Chandra grant GO8-9072X awarded 
through the Harvard College Observatory. S.B. is sup- 
ported in part by a CIFAR Junior Fellowship. C.O.H. 
and I.H.S. are supported by NSERC. This research has 
made use of the NASA Astrophysics Data System (ADS) 
and software provided by the Chandra X-ray Center 
(CXC) in the application package CI AO. 

Facilities: CXO 



REFERENCES 



Abdo, A. A., et al. 2009, Science, 325, 848 

Archibald, A. M., et al. 2009, Science, 324, 1411 

Archibald, A. M., Kaspi, V. M., Bogdanov, S., Hessels, J. W. T., 

Stairs, I. H., Ransom, S. M., McLaughlin, M. A., & Lorimer, 

D. 2010, ApJ, 722, 88 
Arons, J. & Tavani M. 1993, ApJ, 403, 249 
Bassa, C. G., Pooley, D., Homer, L., Verbunt, F., Gaensler, 

B. M., Lewin, W. H. G., Anderson, S. F., Margon, B., Kaspi, 

V. M., & van der Klis, M. 2004, ApJ, 609, 755 
Becker, W. & Aschenbach, B. 2002, Proceedings of the 270. 

WE-Heraeus Seminar on Neutron Stars, Pulsars, and 

Supernova Remnants, Eds. W. Becker, H. Lech, J. Triimper, 

p. 64 

Becker, W., Swartz, D. A., Pavlov, G. G., Eisner, R. F., Grindlay, 
J. E., Mignani, R., Tennant, A. F., Backer, D., Pulone, L., 
Testa, v., & Weisskopf, M. C. 2003, ApJ, 594, 798 

Begin, S., Thesis submitted to the Faculty of Physics, University 
of British Columbia, 2006 

Bogdanov, S., Grindlay, J. E., & van den Berg, M. 2005, ApJ, 
630, 1029 

Bogdanov, S., Grindlay, J. E., Heinke, C. O., Camilo, F., Freire, 

P. C. C, & Becker, W. 2006, ApJ, 646, 1104 
Bogdanov, S., Rybicki, G. B., & Grindlay, J. E. 2007, ApJ, 670, 

668 

Bogdanov, S., Grindlay, J. E., & Rybicki, G. B. 2008, ApJ, 689, 
407 

Bogdanov, S. & Grindlay, J. E. 2009, ApJ, 703, 1557 
Bogdanov, S., van den Berg, M., Heinke, C. O., Cohn, H. N., 

Lugger, P. M., & Grindlay, J. E. 2010, ApJ, 709, 241 
Broos, P. S., Townsley, L. K., Feigelson, E. D., Getman, K. V., 

Bauer, F. E., & Garmire, G. P. 2010, arXiv: 1003.2397 
Cameron, P. B., Rutledge, R. E., Camilo, F ., Bildsten, L., 

Ransom, S. M., & Kulkarni, S. R. 2007, ApJ, 660, 587 
Camilo, F., Lorimer, D. R., Freire, P., Lyne, A. G., & 

Manchester, R. N. 2000, ApJ, 535, 975 
Camilo. F. & Rasio, F. A. 2005, ASP Conf. Ser. Vol. 328; Binary 

Radio Pulsars, eds. F. A. Rasio & I. H. Stairs (San Francisco: 

ASP), p. 147 

Campana, S., D'Avanzo. P., Casares, J., Covino, S., Israel, G. L., 
Marconi, G., Hynes, R., Charles, P., Stella, L. 2004, ApJ, 614, 
L49 

Cash, W. 1979, ApJ, 228, 939 
Cusumano, G., et al. 2003, A&A, 410, L9 

D'Amico, N., Lyne, A. G., Manchester, R. N., Possenti, A., & 

Camilo, F. 2001, ApJ, 548, L171 
D'Amico, N., Possenti, A., Fici, L., Manchester, R. N., Lyne, A. 

G., Camilo, F., & Sarkissian, J. 2002, ApJ, 570, L89 
Damiani, F., Maggio, A., Micela, G., & Sciortino, S. 1997a, ApJ, 

483, 350 

Damiani, F., Maggio, A., Micela, G., k. Sciortino, S. 1997b, ApJ, 
483, 370 

Davidge, T. J. 1996, ApJ, 468, 641 
Davis, J. E. 2001, ApJ, 562, 575 

Edmonds, P. D., GiUiland, R. L., Camilo, F., Heinke, C. O., & 

Grindlay, J. E. 2002, ApJ, 579, 741 
Eisner, R. F., et al. 2008, ApJ, 687, 1019 



Freire, P. C, Camilo, F., Kramer, M., Lorimer, D. R., Lyne, 
A. G., Manchester, R. N., & D' Amico, N. 2003, MNRAS, 340, 
1359 

Freire, P. C. C, Ransom, S. M., Begin, S., Stairs, L H., Hessels, 
J. W. T., Frey, L. H., Camilo, F. 2008, ApJ, 675, 670 

Freire, P. C. C, Wolszczan, A., van den Berg, M., Hessels, 
J. W. T. 2008, ApJ, 679, 1433 

Fruchter, A. S., Stinebring, D. R., & Taylor, J. H. 1988, Nature, 
333, 237 

Grindlay, J. E., Camilo, F., Heinke, C. O., Edmonds, P. D., 

Cohn, H., & Lugger, P. 2002, ApJ, 581, 470 
Harding, A. K. & Muslimov, A. G. 2002, ApJ, 568, 862 
Harding, A. K., Usov, V. V., Muslimov, A. G. 2005, ApJ, 622, 531 
Harris, W. E. 1996, AJ, 112, 1487 

Heinke, C. O., Wijnands, R., Cohn, H. N., Lugger, P. M., 

Grindlay, J. E., Pooley, D., & Lewin, W. H. G. 2006, ApJ, 651, 
1098 

Heinke, C. O., Jonker, P. G., Wijnands, R., Taam, R. E. 2007, 
ApJ, 660, 1424 

Heinke, C. O., Jonker, P. G., Wijnands, R., Deloye, C. J., & 

Taam, R. E. 2009, ApJ, 691, 1035 
Hong, J., van den Berg, M, Schlegel, E. M., Grindlay, J. E., 

Koenig, X., Laycock, S., & Zhao, P. 2005, ApJ, 635, 907 
Huang, R. H. H., & Becker, W. 2010, A&A, 510, 67 
Kargaltsev, O., Pavlov, G. G., & Garmire, G. P. 2006, ApJ, 636, 

406 

Kargaltsev, O., & Pavlov, G. G. 2008, AIPC, 983, 171 

Kirsch, M. G. G., et al. 2006, A&A, 453, 173 

Lyne, A. G., Brinklow, A., Middleditch, J., Kulkarni, S. R., 

Backer, D. C. 1987, Nature, 328, 399 
McCUntock, J. E., Narayan, R., & Rybicki, G. B. 2004, ApJ, 615, 

402 

Pallanca, C, Dalessandro, E., Ferraro, F. R., Lanzoni, B., Rood, 

R. T., Possenti, A., D'Amico, N., Freire, P. C, Stairs, L, 

Ransom, S. M., Bgin, S 2010, ApJ, 725, 1165 
Paltani, S. 2004, A&A, 420, 789 
Pavlov, G. G. & ZavUn, V. E. 1997, ApJ, 490, L91 
Pavlov, G. G., Kargaltsev, O., Garmire, G. P., &; Wolszczan, A. 

2007, ApJ, 664, 1072 
PeUizzoni, A., et al. 2009, ApJ, 695, L115 
Predehl, P. & Schmitt, J. H. M. M. 1995, A&A, 293, 889 
Ransom, S. M., Hessels, J. W. T., Stairs, I. H., Freire, P. C. C, 

Camilo, F., Kaspi, V. M., Kaplan, D. L. 2005, Science, 307, 892 
Ray, P. S., Wolff, M. T., Demorest, P., Cognard, I., Backer, D. C, 

Wood, K. S. 2008, AIPC, 983, 157 
Rutledge, R. E., Fox, D. W., Kulkarni, S. R., Jacoby, B. A., 

Cognard, I., Backer, D. C, & Murray, S. S. 2004, ApJ, 613, 522 
Shawl, S. J. & White, R. E. 1986, AJ, 91, 312 
Smith, R. K., Edgar, R. J., & Shafer, R. A. 2002, ApJ, 581, 562 
Stappers, B. W., Gaensler, B. M., Kaspi, V. M., van der Klis, M., 

& Lewin, W. H. G. 2003, Science, 299, 1372 
Venter, C, Harding, A. K., & Guillemot, L. 2010, ApJ, 707, 800 
Webb, N. A., Olive, J. -F., & Barret, D. 2004, A&A, 417, 181 
Zavlin, V. E. & Pavlov, G. G. 1998, A&A, 329, 583 
ZavUn, V. E. 2006, ApJ, 638, 951 
ZavUn, V. E. 2007, Ap&SS, 308, 297 



